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ABSTRACT 

Context. The CoRoT mission is in its third year of observation and the data from the second long run in the galactic centre direction 
are being analysed. The solar-like oscillating stars that have been observed up to now have given some interesting results, specially 
concerning the amplitudes that are lower than predicted. We present here the results from the analysis of the star HD 170987. 
Aims. The goal of this research work is to characterise the global parameters of HD 170987. We look for global seismic parameters 
such as the mean large separation, maximum amplitude of the modes, and surface rotation because the signal-to-noise ratio in the 
observations do not allow us to measure individual modes. We also want to retrieve the stellar parameters of the star and its chemical 
composition. 

Methods. We have studied the chemical composition of the star using ground-based observations performed with the NARVAL 
spectrograph. We have used several methods to calculate the global parameters from the acoustic oscillations based on CoRoT data. 
The light curve of the star has been interpolated using inpainting algorithms to reduce the effect of data gaps. 

Results. We find power excess related to p modes in the range [400 - 1200] //Hz with a mean large separation of 55.2 + 0.8 //Hz with 
a probability above 95 % that increases to 55.9 + 0.2 //Hz in a higher frequency range [500 - 1250] //Hz and a rejection level of 1 %. 
A hint of the variation of this quantity with frequency is also found. The rotation period of the star is estimated to be around 4.3 
days with an inclination axis of i = 50° ^j. We measure a bolometric amplitude per radial mode in a range [2.4 - 2.9] ppm around 
1000 //Hz. Finally, using a grid of models, we estimate the stellar mass, M = 1.43 + 0.05 M Q , the radius, R = 1.96 + 0.046 R Q , and the 
age ~2.4 Gyr. 

Key words. Asteroseismology - Methods: data analysis - Stars: oscillations - Stars: individual: HD 170987 



1. Introduction 

During the present decade the number of confirmed solar- 
like pulsators - those with acoustic modes excit ed by turbu- 
lent motions in the near-surface convection (e.g. Christensen-| 



Dalsgaard 2004 and reference therein) - has increased enor- 
mously thanks, first, to th e growing number of ground-based 
observing campaigns (e.g. Bedding & Kjeldsen 2 007 [ |Aren toft 
|et al.|2008] ), and second, to the high-precision photometry mea- 
surements provided by space instr umentation such a s WIRE 
(Wide-Field Infrared Explorer, e.g. |Bruntt et al.]|2007]), MOST 
(Microvariab ility and Oscillation s of STars, Walker et al.|2003| ) 



and CoRoT (Baglin et al. 2006 ). The latter has been providing 
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* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil, ESA (RSSD and Science Programme), 
Germany and Spain. 



data with an unprecedented quality both in terms of photomet- 
ric precision and in terms of uninterrupted observation lengths. 
CoRoT h as already observed several main-sequence solar-like 
pulsators (Mic hel et al.||2008] ) while it has enabled to resolve 
the individu al mode s of the oscillations spectra of several F stars 
([Appourchaux et al.|2008||Barban et al.|2009[|Garcia et al.|2009| 
Mosser et al.|2009b| ) and G stars ( [Deheuvels et al.|2010[ Ballot 
et al. in prep.). At least, it allowed to derive a large spacing for 
the faintest targets ( [Garcia et al.|2009[|Mosser et al.|2009b] ). The 
measurement of these seismic parameters already offer a valu- 
able tool for accura te determinatio ns of radii (e.g. |Stello et aL 
2009 ) and ages (e.g. Creevey 2008) of stars, which are specially 



interesting for better understand stellar evolution as well as to 
characterise stars hosting planets. All of these CoRoT obser- 
vations ar e the starting po i nt for a better unders tanding of the 
structure ( [Piau et al.||2009[ |Deheuvels & Michel||2009|) and the 
surface dynamics ( Ludwig et al. 2009; Samad fet al.|2009b"la| of 
this class of stars. 
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The Kepler mission, successfully launched in March 7, 2009 
(Borucki et al. 2009), will also contribute to this field by ob- 
serving stars on very long runs (4 years). The quality of the 



first data on stars showing solar-like oscillations ([B edding et al. 
2010a[ piaplin et al.pilOl |Hekker et alT]|2010b[ jStello et al. 



20 10) promises the asteroseismology a bright f uture on the 



study of stellar interiors and dynamical processes ( Christensen- 
|Dalsgaard & Houdek|2009l|Samadi|20"09l|Mathis|2009| ). 

In this paper we present results about a star recently observed 
by CoRoT, HD 170987 (or HIP 90851). This target is a well- 
known double star where c omponents are separated by 0.7" (e.g. 
Dommanget & Nys 2002 ). The main star is a F5 dwarf star with 
a magnitude my ranging from 7.4 to 7.7 in the literature, while 
the second component has a magnitude around 8.5. 

We start by reporting the latest spectroscopic results ob- 
served by the NARVAL spectrograph (in S ect. 2), which shows 
that this star is very similar to Procyon (e.g.[Allende Pr ieto et aT] 
2002). Then we describe the observations done with CoRoT dur- 
ing 149 days and the interpolation done in the data gaps of the 
light curve (Sect. 3). In Sect. 4 we infer the surface-rotation pe- 
riod of the star from the detailed analysis of the low-frequency 
region of the power spectrum and then, we obtain the global 
properties of the acoustic modes and of the star, respectively in 
Sect. 5 and Sect. 6. We finish in Sect. 7 with a discussion of the 
results and in Sect. 8 with the conclusions of the paper. 



2. Spectroscopic results 

We have observed HD 170987 using the NARVAL spectrograph 
on the 2-m class Bernard Lyot Telescope at the Pic du Midi 
Observatory. We acquired one spectrum on each night of 2009 
July 7, 8 10, 11 and 13. The spectra were co-added to obtain a 
signal-to-noise ratio in the continuum of S /N « 900 in the range 
5800-6400 A. As an example of the quality of the spectrum, a 
small section is shown in Fig.[T] 

We analyse d the spectrum using the semi-auto matic software 
package VWA ( [Bruntt et al.|2004[ |Bruntt|2009] ) wh ich uses at- 
mospheric m odels interpolated in the MARCS gri d ([Gustafsson 
et al.|2008 ) and atomic parameters from VALD (Ku pka et al. 
1999b. The abundances for 239 lines were calculated itera- 



tively by fitti ng synthetic profiles to th e observed spectrum us- 
ing SYNTH (Valenti & Piskunov|1996] ). For each line the abun- 
dances were calculated differentially with respect to the same 
line in a solar spectrum from |Hinkle et al.] ( [20 00). The atmo- 
spheric parameters (T e s, logg, [Fe/H], and Vmicro) were deter- 
mined by requiring that Fe lines gave the same abundance in- 
dependent of equivalent width (EW), excitation potential (EP) 
or ionization stage. Only weak lines were used (EW < 90 mA) 
for this part of the analysis as they are sensitive to adjust the 
parameters, while stronger lines (EW < 140 mA) were used for 
the calculation of the final mean abundances. The uncertainties 
on the atmospheric parameters and the abundances were deter- 
mined by perturbing the best-fitting atmospheric parameters, as 
described by |Bruntt et"aL] ( [2008] ) . 



For F-type stars like HD 170789, non-LTE (local thermody- 
namical equilibrium) effects cannot be ignored. Our model of 
atmospheres assumes LTE and we therefore systematically un- 
derestimate the abundance of Fe i while Fe n is nearly unaffected. 
We followed the approach of |Bruntt et al.|(|2008|) and corrected 
the Fei abundance using the results of |Ren tzsch-H olm ( |1996| ). 
The effect is [Fei/H] NLTE = [Fei/H] LTE + 0.046 dex. This af- 
fects our estimate of the su rface gravity logg, as was discussed 
by Fuhrmann et al. (1997]) in the case of Procyon, which has 
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Fig. 2. Two synthetic profiles with different surface gravity are 
compared to the observed (black) spectrum (logg = 4.35 (green 
line) and 3.95 dex (red line)). The hatched regions are used for 
renormalizing the spectrum. The solid rectangles mark the re- 
gions where the x 2 is computed to determine the optimal value 
of log g. 



a very similar spectral type. Following their approach we used 
the Mgib lines at 5173 and 5184 A to determine logg. We first 
adjusted the Van der Waals broadening constants from VALD 
in order to fit the solar spectrum for the canonical value of 
logg = 4.437 dex. We then determined the Mg abundance in 
HD 170987 from two weak Mg lines and fitted the wings of 
the two strong Mgib lines. An example of the fit is shown in 
Fig. [2] where the observed spectrum is compared to two syn- 
thetic spectra with logg = 4.35 and 3.95 dex. The hatched re- 
gions near the border of the plot are used to normalize the spec- 
trum and the grey- shaded rectangles show the sections of the 
spectrum used to compute the x 2 value needed to determine the 
best value of log g. The mean value of the two Mg i b lines gives 
logg = 4.35 ± 0.22 dex. This value is slightly higher but in 
agreement with the more precise value determined from the anal- 
ysis of Fe i + ii lines: log g = 4.20 ± 0.05 dex. Although the latter 
analysis is affected by our adopted correction for the NLTE ef- 
fect, we use it as our final estimate of the surface gravity. 

From our analysis of the spectrum of HD 170987 we deter- 
mine the following parameters: - 6540 ± 36 + 70 K, log g = 
4.20 ±0.05 ±0.10 dex, and v micro = 1.70 ± 0.10 ± 0.05 km s" 1 . 
We give two uncertainties for these values, the first being the 
intrinsic error and the second being the estimated systematic er- 
ror. The systematic errors were estimated from a large sample of 
F5-K1 type stars (Bruntt et al., in prep.). These two uncertainties 
must be added quadratically (e.g. T q q = 6540 ± 80 K). For the 
mean metallicity we use the Fen lines to get [M/H]=-0.15 ± 
0.06 dex. In Table Q] we list the mean abundances of 16 ele- 
ments compared to the Sun ( [Grevesse et al. [2007 ), which are 
also shown in Fig. [3] Finally, for the projected rotational veloc- 
ity we find v sin/ = 19.0 ± 1.5 km s" 1 , which is significantly 
greater than for Procyon. 

We have also generated a mean profile from the spectrum 
using the Least-Sq uares Deconvolution Technique developed by 
Donati et aL] ( |1997| ). This shows evidence that the spectral lines 
are slightly asymmetric, indicating that the spectrum is contam- 
inated by a fainter star. The effect on the determined spectro- 
scopic parameters are only weakly affected by this. 
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Fig. 1. Comparison of the observed (black) and computed (blue) spectrum of HD 170987. The four solid vertical lines mark the 
lines used in the spectral analysis while the vertical dotted lines show other well-known lines. 
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Fig. 3. Abundance pattern of HD 170987 for 16 elements. Circles and box symbols are used for the mean abundance from neutral 
and singly ionised lines, respectively. Notice the high abundance of lithium. The yellow horizontal bar marks the mean metalicity 
with 1-cr uncertainty range. 



Table 1. Chemical composition of HD 170987 measured relative 
to the Sun (AA). n is the number of spectra lines used. 
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To get a first estimate of the seismic parameters of this 
star, we can use the value for the effective temperature found. 
The Hipparcos parallax combined with a bolometric co rrection 
of -0.01 ±0.01 derived using T e s in the calibration of |Flower" 
(1996) gives L/L© = 6.41 ±1.2 (neglecting interstellar redden- 
ing) and consequently R/Rq = 2.1 ±0.2. A rough compari- 
son with solar-metallicity evolutionary tracks gives a mass of 
M/Mq = 1.4 ± 0.2, which finally yields an estimate of the posi- 
tion of maximum amplitude, v max = 910 + 50 //Hz and the mean 
large separation, Ay ~ 52.7 ± 4//Hz (Kjeldsen & Bedding 1995 ). 
Compared to the other publis hed CoRoT F stars (e.g. |Barban 
et al.|2009| [Garcia et al.|2009)), HD 170987 is the refore notably 



more similar to Procyon (e.g. |Arentoft et al.|20 08 and references 
therein). 



In the catalog of the Geneva-Copenhagen Survey of Solar 
neighbourhood III (Holmb erg et al.|2007] ), based on Stromgren 
photometry, a valu e of my =7. 46 was given . The new reduction 
of Hipparcos data ( [van Lee uwen 2007a |b) gives a parallax of 
ti- 11.21 ± 1.02 mas, corresponding to a distance of 89 ± 9 pc, 
leading to an absolute magnitude of My- 2.71. 



3. CoRoT photometric observations 

HD 170987 has been observed by CoRoT in the seismic field 
during 149 days starting on April 11, 2008 till September, 7, 
2008, during the second long run in the galactic centre direc- 
tion (LRc02). In Fig. [4] we show the N2 light curve (in ppm) 
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corrected for the ageing of the CCD by dividing by a third-order 
polynomial fit. We use time series that are regularly spaced in the 
heliocentric frame (i.e. the so-called Helreg level 2 data, Samadi 
et al. 2007a) with a 32 s sampling rate. The overall duty cycle of 
the time series is 89.7%. 



2000 



1000 



E 
a 

3 



-1000 



-2000 




50 100 
Time (days) 



150 



Fig. 4. N2-Helreg time series of the CoRoT target HD 170987 
corrected from the ageing of the CCD. For plotting purposes 
only every fifth point is shown. 



To compute the power spectrum density (PSD) we used a 
standard fast Fourier transform algorithm and we normalized it 
as the so-called one-sided power spectral density ( Press et al. 
[1992]). The resulting PSD is shown in Fig. [5] 

The existing gaps are mostly due to the crossing of the South 
Atlan tic Anomaly (S AA) (see for more details lAuvergne et al. 
2009). The influence of the perturbations of SAA crossings is 



not completely attenuated by the standard correction procedure 
and thus yields alias peaks at multiples of the CoRoT orbital fre- 
quency of 161.7 //Hz (vertical dotted lines in Fig. [5] (top)). For 
each of them a daily modulation can be observed at ± 1 1 .57 jjHz 
and twice this frequency (see the inset of Fig. [5] (top)). In par- 
ticular, the region where we expect to find the acoustic-mode 
power excess (centered around 900 jjHz) is completely domi- 
nated by this peak structure. Thus, we need to interpolate the 
data in the best way to reduce these effects as much as po ssible as 
it was already done in HD 1757 26 (|Mosser et al.|20Q9bj). For that 
star the algorithm proposed by ( Holdschne ider et al.|1989| ) was 
used giving excellent results. In our case, we adopted a different 
strategy based on the "inpainting" technique using a Multi-Scale 
Discrete CosineTransform (see Appendix A for a more detailed 
explanation on the algorithm). This algorithm has been t ested 
and the results have been verified (see for further details Sato 
et al.|2010 ). Moreover, different groups have also used other in- 



terpolation techniques (linear, based on wavelets, etc) leading to 
results for the global parameters of the modes (such as the mean 
large separation) that agree within the error bars. 



4. Surface rotation 

An independent measurement of the surface rotation can be ob- 
tained through the signature left by the active regions crossing 
the visible stellar disk. It provides a first estimation of the stel- 
lar rotation. Moreover, combined with the measured v sin i = 
19.0 ±1.5 km s _1 (see Sect. 2) and assuming a given radius of the 




100 
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Fig. 5. PSD (in units of ppm 2 /yuHz) of the raw N2-Helreg time 
series (top) and after interpolating the data gaps by a Multi-Scale 
Discrete Cosine Transform as explained in the text (bottom). The 
vertical dotted lines represent the frequency of the CoRoT orbit 
at 161.7 //Hz and the first 10 harmonics. The inserts are a zoom 
around the second orbital harmonic, normalized by its standard 
deviation, showing daily modulation with peaks at ±11.57 /jHz 
and twice this frequency. 



star, we can also estimate the inclination axis of the star, which 
can be very helpful to better co nstrain the fitting techniques to 
characterise the acoustic modes ( Gizon & Solanki 2003 ; Ballot 
|etal.|2006||2008] ). 

The light curve of HD 170987 shows a periodic modulation 
of about 4 days (see Fig. [4]), which could be due to the combina- 
tion of the stellar magnetic activity (star spots) and the rotation 
of the star. It is interesting to analyse this periodicity in a more 
detailed way by studying the low-frequency range of its power 
spectrum. This is shown in Fig. [6] where a clear peak appears at 
2.70 //Hz. This frequency corresponds to a period of 4.3 days, 
which can be associated to the rotation period of the surface of 
the star. 

At lower frequency -around 0.32 //Hz- another structure ap- 
pears which could be related to low-frequency noise. Indeed, 
considering that the radius is estimated to ~ 2 R Q and that the 
observed v sin / is of ~ 19kms _1 , the second peak is very un- 
likely to be produced by the rotation of the stellar surface. 

Spot modelling of the light curve wa s perfor med according 
to the method exposed in Mos ser et al.[ (2009a). Since all the 
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between the wavelet with a given frequency (y-axis) along time 
(x-axis). As showed in Math ur et al.| ([2010), the strength of this 
method is to see the evolution with time of the power, as well 
as to resolve the uncertainty between the fundamental period 
and the first harmonic that could happen, as in the solar case 
(Math uret al.|20 08). The left panel of the plot shows the pres- 
ence of power in the wavelet spectrum for periods between 4 
and 5 days. This power excess is observed all along the 149 days 
of observation. When collapsed over the observation time (right 
panel) we see the accumulation of power at 4.3 days well above 
the 95% confidence level limit. 
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Fig. 6. Zoom on the low-frequency part of the PSD. The spec- 
trum has been computed with a classical FFT using the full time^ 
series. ^ 

stars analyzed in that paper, with a spectral type close to HD <u 
170987, show evidence of differential rotation, we used the fit- §r 



ting method D (Table 2 of |Mosser et al. 2009a| ) with a fixed value fc 
of the differential rotation profile. This method has proven to be 
efficient for the determination of the rotation period and of the 
spot lifetime. Measurements of the star inclination and differen- 
tial rotation are much less precise. 



for a 60-day long 
ling is not able to 



The best fit model is presented in Fig 
sample. Like similar type of stars, the mode! 
reproduce the sharpest features of the light curve. Fitting them 
would require the introduction of too many spots, which lowers 
the global precision of the parameters of the best fit. This fit is 
obtained for an inclination of 40^°. The equatorial rotation pe- 
riod is 4.30 ± 0.05 days and the differential rotation rate 5^° %, 
for a mean rotational period close to 4.40 d. The spot lifetime, 
2.5^-^ days is significantly shorter than one rotation period of 
the star. 
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Fig. 7. Spot modelling of HD 170987. The dots in the upper 
panel represent the data binned every CoRoT orbit (6184 s) and 
the solid curve is the best fit model. Dots in the bottom panel 
represent residuals. Vertical grey lines indicate the mean rota- 
tion period. 
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Fig. 8. Left panel: Wavelet power spectrum for HD 170987 as a 
function of the frequency of the Morlet wavelet (y-axis) and time 
(x-axis). The black grid represents the cone of influence. Right 
panel: Global power spectrum as a function of the frequency of 
the wavelet. The dotted line is the confidence level correspond- 
ing to a 95% probability. 



The reliability of the result depends on two factors: the cone 
of influence, which is related to the fact that the periodicity has 
to be at least one quarter of the total length of the time- series and 
a confidence level corresponding to a 95% of probability that the 
peak observed is not due to noise. Thus, using this methodology 
we determine that the rotation period of HD 170987 is -4.3 days 
(or 2.7 //Hz), which is compatible with the analysis of the low- 
frequency PSD described above. 

It is very clear that the power present around 2.7 /jHz is con- 
centrated between the 30th and the 70th day (black shaded re- 
gion) as well as during another small period around the 100th 
day. This agrees with the results obtained with the PSD. 

Given the value of the rotation period and the v sin i, we can 
deduce the inclination angle i = 50° 



+20 
-13* 



5. Global parameters of the acoustic modes 

Using different methods, we have estimated the global parame- 
ters of the acoustic modes of HD 170987, such as the location 
of the p-mode bump, the mean large spacing, and the maximum 
bolometric amplitude per radial mode. 



We have used the wavelet te chnique to confirm the rot ation 
period of the surface of the star ( [Torrence & CompofT9 98 ). We 
have used the Morlet wavelet (a moving Gaussian envelope with 
a varying width) to produce the wavelet power spectrum shown 
in Fig. [8] This spectrum is a representation of the correlation 



5.1. P-mode excess power in the PSD 

To evaluate the possible presence of p-mode oscillations in the 
spectrum we start by analyzing a smoothed version of the PSD. 
We smoothed the PSD in the range 200 to 8000 //Hz with a 
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Gaussian running mean with a width of 200 //Hz. The reason 
not to include the low-frequency part of the spectrum is that this 
region of the PSD mainly bears the signatures of the slow trends 
in the data and of stellar activity (rotation), which is already ex- 
plained in Sect. [4] The smoothed PSD is shown in Fig. [9] Two 
clear bumps are seen in the spectrum, one around 300 and one 
around 1000 //Hz. The bump around 1000 jjHz could originate 
from the p-mode oscillations (as it corresponds to the expected 
Vmax) but it could also be a signal related to granulation. Up 
to now, the background of all the precedent CoRoT solar-like 
targets (HD 49933, HD 175726, HD 181420, HD 181906 and 
HD 49385) has been correctl y characteriz ed by a single compo- 
nent in the Harvey model (Harvey 1985) that has been associ- 
ated to granulation. However, when comparing the spectrum of 
HD 170987 to the one of the Su n obtained using the VIRGO in- 
strument (Frohl ich et al.|1995] ) and to solar-type stars observed 
by the Kepler satellite (Chaplin et al. 2010) we see that in these 
cases a second Harvey power law has been used to take into ac- 
count some extra power located around 1000 jjHz (time-scale 
close to 80 s). Based on the smoothed PSD we can therefore not 
conclude that the spectrum shows p-mode oscillations while we 
see some excess power centered at ~ 1000 jjHz. 

For the precise analysis of the acoustic background, different 
strategies have been followed to study the bump at 1000 //Hz. We 
start by assuming that the bump around 300 jjHz originates from 
granulation and the one at ~ 1000yuHz could originate from a 
shorter convective scale or by faculae. 

The faculaes are the bright points seen on visual solar im- 
ages often close to the dark sunspots. They are due to changes in 
the opacity caused by a strong magnetic field which means that 
we see the inside rather than the surface of the granulation and 
as the temperature is higher inside the granulation cells than at 
their surfaces the facul aes appear brighter than there surround- 
ings ( [Keller et al.|2004| ). The carateristic timescale of the faculae 
is shorter than the time-scale of the granulation cells because the 
faculaes only sample the edges of the granulation cells. 

Following Harvey (1985) we model the smoothed PSD in 
the range [200 - 8000] //Hz with a background model contain- 
ing two components plus the photon-noise contribution (c) that 
dominate the PSD at high frequency, i.e. 



PSD(v) = c + 



Act 2 t 

^ u gran * gran 
1 + (27TT grm v) a 



4<T facu Tfacu 
1 + (27TT facu v)^c 



(1) 



Here, <x gran and <Tf aC u are the amplitudes of the background sig- 
nal of granulation and the second convective component or fac- 
ulae respectively, r is the time-scale, and a gran and <2f acu are two 
constants. Because we are modelling a smoothed version of the 
PSD we can assume that the error between the model and the 
observations are normally distributed. We can therefore fit the 
model to the observations by means of least squares. Using the 
robust non-linear least squares curve fitting IDL package MPFIT 
(http://www.physics.wisc.edu/~craigm/idl/idl.html) we obtain a 
granulation time-scale of 383 + 28 s and a second component 
with a time-scale of 1 13 ± 51 s. 

We then assumed that the second bump, around 1000 yuHz, is 
caused by p modes. Thus, for the background fitting, we use the 
first two terms of equation (2) - the photon noise and one Harvey 
model component- and fitted a Gaussian function to the second 
bump. This led to a very similar fit compared to the previous 
one, with the same parameters for the granulation time-scale (see 
Fig. [9]). Consequently, both assumptions could be correct and we 
cannot disentangle which one is the best one at this stage of the 
analysis. 




500 1000 5000 

Frequency [/^Hz] 

Fig. 9. Smoothed PSD (continuous solid line) overlaid by the 
background fitted model assuming two Harvey-model compo- 
nents (dotted lines) and a single Harvey model plus a p-mode 
excess (dashed lines). The individual components of the fitted 
background are also shown. The PSD have been smoothed with 
a Gaussian running mean with a width of 200 jjHz. 



5.2. Estimating the mean large separation 

In order to determine whether HD 170987 presents solar-like 
oscillations, we have applied several methods. One of them con- 
sists of searching for the signature of the mean large separa- 
tion of a solar-like oscillating sign al in the autocorre lation of 
the time series. As proposed by Roxburgh & Vorontsov (2006), 
in a refined application of the Wiener- Khinchine theorem, this 
autocorrelation is calcula ted as the Fourier spectrum of the win- 
dowed Fourier spectrum. |Roxbu rgh (2009) subsequently used a 
narrower window to look for the variation of the large separation 
with frequency. 

Mosser & Appourchaux ('2009') have shown how to optimize 
the method and to determine its reliability. They have scaled the 
autocorrelation function according to the noise contribution, in 
order to statistically test its significance. Then, when the enve- 
lope autocorrelation function (EACF) gives a signal above a de- 
fined threshold level, the null hypothesis can be rejected, and a 
reliable large separation can be derived. For a blind analysis of 
the mean value, (Ay), of the large separation and a rejection of 
the null hypothesis at the level 1 %, the threshold level is fixed 
to 8. This mean value is determined in a frequency range cen- 
tered on the frequency vracf where the EACF reaches its maxi- 
mum amplitude, with a filter of width equals to the full width at 
half-maximum of the mode envelope. Note that veacf is close to 
Vmax, where the mode envelope reaches its maximum amplitude, 
but they are not equal. The difference between veacf and v max is 
due to the varying mode lifetime. 

The method has shown to be efficient in low signal-to-noise 
cases, such as the CoRoT target HD 175726 (Mo sser et al.| 
2009b]) or in the Kl V target with a he ight-to-background ratio <R 
as low as 2.5 % ( |Gaulme et al.|2009| ). The ratio % measures the 
maximum smoothed height of the modes compared to the back- 
ground. HD 170987 is also a challenging target, with % ^ 4 %. 
For HD 170987, with an optimized filter of width 580yuHz, the 
EACF peaks at the frequency 930yuHz, at a value of 19.2, and 
gives a mean large separation of 55.9 ± 0.2yuHz. 
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Fig. 10. Maximum relative power of the highest peak in the PS2 
around half of the large separation, as a function of the cen- 
tral frequency of the sliding box taken in the PSD. The dashed 
line represents the threshold corresponding to a 95% confidence 
level. The horizontal error bars represent the boxes in which the 
PSD has been sliced (600 jjHz in this case). 



Other methods have been used giving a mean large separa- 
tion b etween 54 and 56 jMz ( [Hekker et al.|2010a|[Mathur et al] 
|20T0] ). 



Indeed, using the method described in Mathu r et al. ( |2010 ), 
we have calculated the power spectrum (PS) of the power spec- 
trum (PS2) between 100 and 10000 //Hz. We select the high- 
est peak and we assume that it corresponds to half of the mean 
large separation. Then, we cut the original PSD in boxes of 
600 //Hz, shifted every 60 /jHz and we compute the PS2 of each 
box looking for the highest peak close to the value found origi- 
nally. In Fig. 10 we present the power of the highest peak found 
in each box (normalised by the cr of the PS 2). The horizontal 
error bars mark the frequency range covered by each box. We 
have obtained that the frequency-range of the p modes is [400 
- 1200] //Hz with <Av> = 55.2 ±0.8 fiUz with 95% confidence 
level. 

Finally, we took into account that Ay depends on frequency 
and we computed the PS2 of a power spectrum with a stretched 
frequency axis (He kker et al. [2010a). The stretching is per- 
formed in such a way as to produce an equidistant pattern of 
peaks on the stretched, as opposed to the original, frequency 
axis. It is done by applying a quadratic correction to the axis val- 
ues. This allows us to measure the linear variation of the Ay with 
frequency (which are in themselves first differences of frequen- 
cies). The PS2 of the stretched power spectrum will therefore 
show a stronger (more prominent) signature of Ay than the PS2 
of the original spectrum. In the PS2 we determine the position of 
the Ay/2 and Ay/4 features, and compute the Bayesian posterior 
probability of the features being due to noise. Using these prob- 
abilities we then compute the posterior weighted centroids of the 
features. The interval with a probability of the feature not being 
due to noise higher than 68.27%, i.e., Icr in a Gaussian distri- 
bution, is used as the uncertainty interval. For this star we could 
only find a Ay/2 feature in the PS2, which was used to determine 
(Ay). 



Let's note that a large separation around 55 fiRz is very close 
to the third harmonic of the orbit (53.9 //Hz) which are still 
present in the inpainted data. Therefore, we have repeated the 
analyses described before over a cleaned version of the power 
spectrum in which we have removed by hand the remaining 
peaks of the orbital harmonics substituting the affected bins by 
a local average of the PSD (see Fig. [TT]). The results remained 



unchanged which means that the pattern is not a consequence of 
an instrumental artifact. 

Using the large separation previou sly computed, we can cal- 
culate the so-called echelle diagram ( Gre c et al. |1983| ). To do 
so, we fold the power spectrum using a large separation of 55. 
//Hz in the region where the p-mode excess has been found [440 
- 1200] //Hz. To reduce the dispersion of points in the spectrum 
we have first smoothed it by a boxcar function of 40 points (3.1 
//Hz). The result is shown in the top, left panel of Fig. [TT] This 



diagram is dominated by the presence of the harmonics of the 
orbit at multiples of 161.7 //Hz, which are seen as an inclined 
line of points starting at 48 //Hz. If we use the cleaned version of 
the PSD (where we have removed most of the peaks related to 
the orbital harmonics) and we do the same analysis (see top right 
panel in Fig. 11), two nearly vertical ridges clearly appear in the 
diagram. Indeed, if we integrate (collapse) the echelle diagram 
in the vertical direction we obtain two well separated ridges (see 
bottom panels in Fig. 11). We can tentatively assign the left-hand 
ridge with even-1 modes and the right-hand one to the even ones 
odd knowing that the signal-to-noise ratio is not enough to go 
beyond that. 

The values of the mean large separation found by the dif- 
ferent groups agree within their error bars and in the following 
sections, we will take (Ay) = 55.2 + 0.8 //Hz for the range [400 
- 1200] //Hz. It is important to note that the error bar given here 
is the internal error of one single method and it does not take into 
account the dispersion in the results from the different pipelines 
and the slightly different ranges used in their computations. 



5.3. Variation of the large separation 

When calculated with a filter narrower than the mode envelope, 
the EACF makes it possible to investigate the variation Av(v) 
of the large separation with frequency. With a filter width equal 
to 5 times the mean large separation, we test the variation at a 
large scale: the large separation increases with frequency at a 
rate dAv/dn of about 0.55 ± 0.15 //Hz per radial order, reaches 
its maximum at 980yuHz and then decreases (Fig.[T2|). 

With a filter width equal to 2 times the mean large separa- 
tion, we can investigate the rapid variation of the large sepa- 
ration. The number of independent measurements allows us to 
show the modulation of the large separation (Fig. [13]). For the 
study of the variation of the large separation with frequency 
Ay(y), the threshold level is lower than the one for the mean large 
separation, benefitting from the fact that the large separation is 
then searched in a pre-determined reduced interval compared to 
the blind analysis. Most of the values of the EACF with such 
a filter are above 4.6, which corresponds to the 1 % rejection 
level for the analysis with a 2- (Ay) wide filter and a search in 
the frequency range allowing 20 % variation with respect to the 
mean large separation. They give error bars lower than 1 yuHz, 
except in the ranges [790 - 900]yuHz and [1180 - 1280] //Hz. 
The EACF values are large enough to enable an unambiguous 
detection of the modulation. The period W of this modulation, 
close to 200yuHz, could be caused primarily by the Hell ion- 
isation zone or it could be related to the internal phase shifts. 
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Fig. 11. Echelle diagrams obtained with the inpainted PSD (left panel) and the cleaned PSD (right panel) as explained in the text. 
They are computed using a folding frequency of 55 /jHz in the region [440 - 1200] /jHz where the excess p-mode power was found. 
Bottom panels are the integrated (collapsed) echelle diagrams along the vertical axis. The dashed line represents the places where 
the harmonics of the orbit are situated. 




0.8 
Frequency 



Fig. 12. Variation of the large separation with frequency. 1-cr 
error bars are indicated. The extent of the Hanning filter used 
by the EACF is indicated in the lower-right corner, the dotted 
line measuring the full- width at half-maximum of the filter, here 
equal to 5 times the mean large separation. The width of the 
white and grey vertical regions is equal to (Ay). 



Following Monteiro & Thompson ( 2005 ), we infer the acoustic 
depth r ^ 1/2W of this region in the range r = 2 300-2 700 s. 



The analysis of the variation of the large separation with fre- 
quency, Av(v), shows a peculiar behavior in the range [790 - 
900]yuHz. One would expect large values of the EACF in this 
range close to v max , which is not the case. Since the variations 
of Av(v) are not greater in this region than in others, one can ex- 
plain the low values of the EACF either by a perturbation caused 
by a mixed mode or by a low value of the mode lifetime. The 
presence of a mixed mode is in fact quite unlikely: such a mode 
with a much longer lifetime than normal p mode should appear 
with a large amplitude, which is not the case. Thus, one has to 
favor the hypothesis of the varying mode lifetime. The compari- 
son of the EAC F with the bolometric amplitude per radial mode 
(see Sect. |5.4| ) indicates a strong increase of the mode lifetime 
below 850yuHz and a strong decrease above 1100//Hz (Fig. [14). 
This behaviour is in agreement with the measurement made in a 
similar star with a very close large separation, such as HD 49385 
( [Deheuvels et al.|2()TU| ), in spite of a different effective tempera- 
ture. However, it is not yet possible to quantify the lifetime. 



5.4. Maximum bolometric amplitude 

To compute the maximum bolometric amplitude per radial mode 
we need to smooth the p-mode hump and correct from the back- 
ground in this frequency range. To do so, we have used differ- 
ent methods (see for an extensive explanation Kjeld sen et al.| 
|2008[ |Hekker et al.|[20T0al |Mathur et al1|2010> . Then to cor- 
rect from the instrumental response function we use the method 
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Fig. 13. Same as Fig. 12 but with a narrow filter, corresponding 
to 2 times the mean large separation. The thin lines correspond to 
frequency ranges within the EACF is belove the threshold level 
for a rejection of the null hypothesis at the 1 % level. The dashed 
lines are the regions perturbed by the presence of orbital har- 
monics. 



( [19951 ) with s = 0.7 ( [Samadi etal.|2007b) . It can be seen 
e in all cases the amplitudes are systematically about 
30 % lower than the theoretical values, the discrepancy appears 
to be significantly larger for HD 170987 (~ 50 %). We note that 
the error bars are larger for HD 170987 due to the rather large 
uncertainty on the parallax (-10 %). In both cases, the discrep- 
ancy is well inside 3-cr. 
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Fig. 14. Comparison of the EACF (dashed line), estimated with a 
5-(Av) large filter, and the bolometric mode amplitude expressed 
in ppm (continuous line). 



developped by Michel 'et al.| (2009). Thus, we have obtained a 
bolometric amplitude per radial mode of Ab i,i=o = 2.7 ± 0.2 ppm 
at -1070 //Hz. However, among the different teams and depend- 
ing on the method and on the way we fit the background, the 
value for the maximum amplitude per radial mode varied within 
the range 2.4 to 2.9 ppm indicating that the error bar of 0.2 ppm 
could be underestimated. So we have estimated a new value for 
the error of 0.6 ppm, which is derived from the scatter of the 
smoothed power spectrum about the background fit outside the 
oscillation range (|Huber et al.|2009|). 



The top panel of Fig. 15 shows the comparison of the am- 
plitude found for HD 170987 with the published values of four 
other C oRoT solar-like targets ( [Michel et al.|2008|[Mosser et al. 
2009b ) as a function of v max . The bottom panel shows the am- 
plitudes compared to theoretical values for each st ar calculated 
using the (L/M) S T~®- 5 scaling relation from Eq. 3 of Kjeldsen & 



Fig. 15. Top panel: Maximum bolometric amplitude per radial 
mode versus v max for published CoRoT F-stars (open symbols), 
the Sun (star symbol) and HD 170987 (filled circle). Bottom 
panel: Maximum bolometric amplitude per radial mode ver- 
sus stellar parameters. Black symbols are the observed values, 
and blue symbols theoretical values using (L/M) 07 . The symbol 
types are the same as in the top panel. 



6. Inferring global stellar parameters 

Using (Av) = 55.2 ± 0.8//Hz, U n = 400yuHz, / max = 
1200 //Hz (see Sect. 5.1 and 5.2), and the spectroscopic infor- 
mation logg = 4.20 ± 0.14 dex, [M/H] = -0.20 ± 0.15 dex 
and r ef f = 6540 ± 80 K (see Sect. [2]), we compare these val- 
ues with stellar models to determine the radius and the mass 
of HD 170987. The stellar models that we use are the Aarhus 
Stellar Evolution Code (ASTEC) coupled with an adiabatic pul- 
sation code (ADIPLS) (Christensen-Dalsgaard 2008a,b). These 
codes need as input the stellar parameters of mass, age, chemi- 
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Table 2. Stellar parameters from modelling 



Parameter 


Value 


Error 


Radius (R ) 


1.960 


0.046 


Mass (M Q ) 


1.43 


0.05 


Age (Gyr) 


2.4 





cal composition, and mixing-length parameter, and return stellar 
observables 5,-, such as radius, effective temperature, and the fre- 
quencies of the oscillation modes. 

The parameters that best describe the observables are ob- 
tained by minimizing &x 2 function; 



X 



M 

Z 

i=\ 



' y_i - Bj 



(2) 



where yi and ^ are the i = 1,2, ...,M observations and errors. 
Here, we have M = 4. The Levenberg-Marquardt algorithm is 
used for the optimization, and this incorporates derivative in- 
formation to guess the next set of parameters that will reduce 
the value of^ 2 . Naturally, an initial guess of the parameters is 
needed and these are obtained from a small grid of stellar evolu- 
tion tracks. 

Because there are few observations and just as many parame- 
ters, there are inherent correlations between mass, age and chem- 
ical composition. To help avoid local minima problems, we min- 
imize the x 2 function beginning at several initial guesses of the 
parameters (mainly varying in mass and age), and these initial 
guesses are estimated from the grids. We therefore obtain sev- 
eral sets of parameters with a corresponding^ 2 value that match 
the observations as best as possible. Choosing the models whose 
X 2 values fall below 3.9 2 (assuming that the errors are Gaussian 
and that we want a 99 % probability), the radius and mass are 
determined by calculating the average value of each, and the un- 
certainty is given by the dispersion in the values divided by 6 (see 
Mat hur et al.|2010| for details). We estimate the age of the star 
by simply matching the "fitted mass" with the closest age from 
the model parameter results, we therefore have no quantitative 
measure of its uncertainty, leading to a very preliminary value. 
The results are given in Table [2] The estimate of the age implies 
a star coming close to the end of the main sequence phase. 

7. Discussion 

Thanks to spectroscopic observations from NARVAL, we have 
been able to estimate the effective temperature as well as the 
chemical abundances of HD 170987 showing that this star has 
a low metallicity (-0.15 ± 0.06 dex) and a high abundance of 
lithium. 

The surface light elements abundances could give us an in- 
teresting clue about the internal rotation of HD 170987. The 
effective temperature of 6540 K sets this star in the so-called 
lithium dip where both 7 Li and 9 Be are depleted (Boesgaa rd & 
King 2002 ). This abundance dip is most likely related to rotation 
and observations show that depletion and r otation are correlated 
( [Boesgaard|[T987] ). |Talon & Charbonnel| ( [2003] ) explained the 
cooler side of the lithium dip (from 6500 to 6200 K) by invoking 
an increasing efficient transport of the internal angular momen- 
tum through gravity waves when the mass decreases from ^1.4 
to « 1.2 M Q . For the more massive stars differential rotation in 
the interior leads to shear instabilities: the radiative envelope is 
mixed deeply and light elements are destroyed through proton 



capture. For the less massive stars there is less differential rota- 
tion and less shear instabilities. Thus no (or far less) depletion 
occurs. 

The lithium dip concerns stars that have an effective tem- 
perature in the range 6400 to 6850 K. At Hyades age this star 
should already have exhib ited a depletion of about 0.5 dex in 7 Li 
([Boesgaard & King 2002). However HD 170987 seems much 
older than Hyades and shows no such depletion: if we combine 
the current solar 7 Li abundance (Asplu nd et al.| [2009) and the 
relative abundance of HD 170987 to the Sun given above we ob- 
tain A( 7 Li) « 2.9 dex. This value is very close to the fraction 
of lithium in the material out of which Population I stars form: 
3.3 to 3.2 dex. Thus the abundance of 7 Li in the atmosphere of 
HD 170987 directly supports the absence of deep mixing and in- 
directly supports the absence of differential rotation with depth. 
At least if some internal differential rotation exists it is not suffi- 
cient to induce deep mixing. 

As confirmed by our analysis, HD 170987 is remarkably 
similar to Procyon, a prime target for previous asteroseismic 
studies using ground-based spectroscopy and space-based pho- 
tometry. For the latter, observations by the MOST satellite re - 
sulted in a null-detection of p-modes (Matthews et al. 2004), 
leading to conclusions that were later critically discussed by 
Bedding et al.| ( [2005 1). These observati ons were followed by 
WIRE photometry ( [Bruntt et al.||2005|) and high er precision 
MOST data collected in 2007 ( [Guenther et al.|2008) , both show- 
ing power excess around the expected p-mode location of lmHz, 
while velocity observations revealed a clear p-mode structure 
(see, e.g., |Arentoft et al.|2008[|Bedding et al.|2Q10b] ). Convection 
models by Guenthe r et al.| ( |2008| ) showed that for Procyon, as 
opposed to the Sun, the granulation timescale can coincide with 
the location of pulsation frequencies. The resulting interaction 
with the acoustic oscillations causing short mode lifetimes were 
identified as possible difficulties for detecting clear peak spac- 
ings in the photometry of Procy on-like stars. Given the similar- 
ity of HD 170987 to Procyon, it is tempting to speculate that a 
similar mechanism is in part responsible for the low visibility of 
the pulsation modes in the CoRoT photometry of HD 170987. 
Alternate mechanisms which could partly explain the low signal 
include the low metallicity (see |Samadi et al.|2009b| ) or the fact 
that some of the light is diluted by a possible binary companion. 



8. Conclusions 

We have analysed the 149 days of the light curve of the CoRoT 
target HD 170987, allowing us to determine its global param- 
eters. We find that the rotation period of the star is around 4.3 
days and the v sin i of 19 km/s, leading to an inclination angle i 
of ~ 50° ^3. By fitting the background, we obtain that the time- 
scale for the granulation is around 383 ±28 s either assuming 
that we have some contribution of faculae or of modes. 

We confirm that power excess between 400 and 1200 //Hz 
is due to acoustic modes, which are characterised by a mean 
large separation (Ay) = 55.2 ±0.8 jjHz with more than 95% 
of probability and a posteriori probability of 68 %. Using the 
EACF method, we find (Ay) = 55.9 ±0.2yuHz in a higher fre- 
quency range [500 - 1250] jjHz with a 1 % rejection. However, 
we recognize that all these confidence levels are dependent on 
our assumption that our spectrum is dominated by white noise 
statistics. It may therefore be that we should be less confident 
as there could be other sources of noise embedded in the sig- 
nal. However, at such a level, the detection seems to be un- 
ambiguous. Moreover, this value is in agreement with the ex- 
pected one derived from the stellar fundamental parameters, 
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which gives 52.7 + 4 //Hz, the large error bar being due to the 
large uncertainty on the parallax. Now, thanks to the detection of 
this pattern, we can be confident that the bump observed around 
1000 //Hz is very likely to be related to acoustic modes. 

We characterise the power excess with a maximum ampli- 
tude of 2.7 ±0.6 ppm at -1070 yuHz, where the uncertainty on 
the amplitude has been obtained by taking into account the scat- 
ter of the smoothed power spectrum about the background fit 
outside the oscillation range. 

Because of the low signal-to-noise ratio, we cannot fit the 
acoustic modes individually. However, from the global seis- 
mic parameters found and the results form the NARVAL spec- 
trohraph, we have estimated that HD 170987 has a mass, 
M = 1.43 ±0.05 M , a radius, R = 1.96 ±0.046 R Q , and an 
age -2.4 Gyr. 

Further studies are needed to have a better comprehension of 
the propagation of waves in stellar atmosphere of stars located 
at this position in the HR diagram. It would also enable us to 
understand why the oscillation-amplitudes measured in intensity 
fluctuations are so small. 
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Appendix A: Inpainting interpolation of the light 
curve 

"Inpainting" refers to methods that interpolate the missing in- 
formation using some priors on the solution. The method that 
we applied to HD 170987 uses a prior of sparsity and was in- 
troduced by Elad et al. (2005). It has already been applied to 
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various fields of astrophysics, such as weak lensing ( Pires et al 
2009) and recently to asteroseismic time series by Sato et al. (in 



prep.). It assumes that there exists a dictionary O - which in our 
case is a Discrete Cosine Transform (DCT) - where the com- 
plete data are sparse and the incomplete data are less sparse. It 
means that there exists a representation a = O r X of the signal X 
in the dictionary <D where most coefficients a\ are close to zero. 
Let X be the ideal complete time series, Y the observed time se- 
ries and M the window function or mask (i.e. M t - 1 where there 
is a valid data point and M t = elsewhere). Thus we have that 
Y = MX. Inpainting consists in recovering X knowing Y and M. 
The solution is obtained by minimizing the following equation: 

min ||O r X||i subject to ^(7 - MX) 2 < a, (A.l) 

i 

where <x stands for the noise standard deviation and where we 
use a pseudo norm, with ||z||i = Yjk I Zk I- To minimize this 
expression, an iterative algorithm was used. 

In this study, we have built the mask M to remove those 
points which are affected by the S AA and also all the other points 
that were already linearly interpolated, according to the status 
flag provided in the N2 files. 

Most of the gaps in the CoRoT time series have a typical time 
scale of several minutes. It is possible to interpolate these typical 
gaps by using an inpainting algorithm based on a DCT decompo- 
sition which is localized in a short time scale. However, there are 
also a few longer gaps which have a time scale of several hours 
in the time series. The inpainting algorithm needs to take into ac- 
count this large variation of the scale of the gaps. To do that, we 
have replaced the simple DCT decomposition by a Multi-Scale 
Discrete Cosine Transform decomposition. With this transform, 
we first decompose the light curve into different scales using a 
wavelet transform. Then, the inpainting of each scale makes pos- 
sible to interpolate the gaps of different time scales. The PSD of 
the inpainted time series is shown in Fig [5] (bottom). The im- 
provement is clearly visible by a reduction of the amplitude of 
the orbital harmonics by more than 90 % in many cases, as well 
as the disappearance of the daily modulation except for the first 
orbital harmonic (see the inset of Fig [5] (bottom) corresponding 
to the second orbital harmonic in which the main peak is reduced 
by 95 % and no signature of the daily modulation is seen). 



